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Abstract 

The effect of tliermal fluctuations on the behaviour of a 4-terminal 
SQUID is investigated. The studied system consists of the Josephson 
4-terminal junction with two terminals short-circuited by a supercon- 
ducting ring and the other two form the transport circuit. The be- 
haviour of a 4-terminal SQUID is controlled by external parameters, 
the applied magnetic flux and the transport current. They determine 
the voltage in the transport channel and the magnetic flux embraced 
in the ring. Within the numerical model the noise-rounded current- 
voltage characteristics and the time dependence of the magnetic flux 
have been calculated. In some region of the control parameters the 4- 
terminal SQUID is in the bistable state with two magnetic flux values. 
The switching of the flux in bistable state produced by the thermal 
noise or by the transport current is studied. 

Keywords: Josephson junction; Multiterminal;SQUID;Thermal fluc- 
tuations 



1 Introduction 



Presently the properties of superconducting microstructures attract theoreti- 
cal and experimental interest. Mostly it is due to the high level of technology 
which permits the fabrication of the ultrasmall Josephson weakly coupled 
systems with well defined geometry. Thus, the well known effects of weak 
superconductivity ^ now can be studied in the situation approaching to 
the idealized theoretical models. From another side, new physics (e.g. un- 
conventional proximity effect, reduced shot noise, Andreev's interferometry) 
appears in small mesoscopic structures f^. The conventional Josephson ef- 
fect in microconstrictions also can be considered as mesoscopic phenomena 
produced by the macroscopic interference in the constriction of the initially 
coherent Cooper pairs wave functions of the two massive superconducting 
banks. The old idea[^ |] of more complicated Josephson weak link struc- 
tures, the so called Josephson multiterminals 0, was realized experimentally 
in ref. |^ . In such structures the quantum interference takes place in the 
intersection of several microbridges connecting superconducting terminals. 
The specific multichannel interference effects were studied theoretically and 
experimentally in the novel superconducting device, the 4-terminal SQUID 
controlled by the transport current [|r^ . In this system we have the interplay 
between the coherent current states in the superconducting ring placed in 
the magnetic field and the Josephson effect in the transport circuit. The 
external transport current / and the applied magnetic flux $e are the con- 
trol parameters which govern the behavior of the 4-terminal SQUID. These 
two independent parameters determine the response of the system: the total 
magnetic flux $ in the ring and the voltage V over the transport channel. It 
was shown |T0| that in the steady state domain (when the transport current 



is less than the critical value /c($e)) the nonlinear parametric coupling leads 
to a bistable state of the SQUID $($e; I) for all values of the self-inductance 
of the ring. Outside the steady state domain (/ > 1^) the resulting flux 
$ becomes an oscillating function of time whose amplitude, frequency and 
shape can be smoothly controlled by the external parameters |1TT| , |12[. 

In the papers |ll|, the theory of a 4-terminal SQUID was developed 
without taking into account thermal fluctuations. The noise in the system 
obviously plays an essential role. In particular it leads to the magnetic flux 
switching in the bistable state, as was directly observed in experiment [p!3| . 

In the present paper we study the effects of thermal fluctuations on the 
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behavior of a 4-terminal SQUID. In Section 2, the basic equations used to 
describe the system are presented. The thermal noise is considered by adding 
the fluctuating voltages to the dynamical equations for the phases of the 
superconducting order parameter. These equations were solved numerically 
following the procedure described in In Section 3 the calculations of 



the noise-affected current-voltage characteristics for several values of applied 
magnetic flux are presented. The time dependence of the magnetic flux in 
the ring stimulated by noise is studied in Section 4. 



2 Equations for the 4-terminal SQUID 

The 4-terminal SQUID which is controlled by the transport current is shown 
in Fig.l. It consists of the Josephson 4-terminal junction with two terminals 
(3 and 4) short-circuited by a superconducting ring and the other two (1 and 
2) form the transport circuit. 

The set of equations describing the behavior of a 4-terminal SQUID in an 
applied magnetic flux and a transport current was obtained in ref.[|lOl- The 



Josephson supercurrent in the jth filament (connecting the centre o with 
the jth bank) expressed in terms of the phases of the superconducting order 
parameter (related to the phase in the center o) has the form 

js ^ > Sm(09,- - Ifiu) (1) 

^ClvJ-c j2 l/Rf^ k=l ^^]^k 
k=l 

Here Aq is the gap in the bulk banks, Rj is the normal resistance of the jth 
filament. Let Vj be the voltage in jth terminal, related to the voltage in the 
centre ( J2 ^jl Rj = 0). In the frame of the heavily damped resistively shunted 

3 

model0 we add to supercurrents (1 ) the normal currents /" = Vj/Rj with 
Vj = 'iPj. Thus for the total current flowing in jth branch we have the 
expression in terms of the phases (pj 

h - AekT. ^ ^^^^ R,R, SmlC^, ^k) + 2e/?, dt + R, ' 

J = 1,2,3,4. 



3 



The positive sign of Ij (2) corresponds to the direction of the current 
from the centre to jth bank. Note that in the resistively shunted model the 
superconducting and normal currents are conserved separately and equations 
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(2 ) satisfy to conservation of the total current, J2 Ij = ^ ■ 

i=i 

The third term in the r.h.s. of equation ( 2) represents the thermal John- 
son noise of the normal resistance Rj. We assume that the time dependent 
voltage noise sources are uncorrelated, each having a white voltage spectral 
density Sy = ^kTRj. 

The values of the currents Ij are determined by the circuit implication of 
the Josephson multiterminal. In the 4-terminal SQUID configuration consid- 
ered here (see Fig.l) they are connected to the external controlling parame- 
ters, the given transport current / and the applied magnetic flux in such 
a way that 

h = -/, h = I,h = ~J, h = J,J= ($e - <^)/L, (3) 

where J is the circulating current in the ring with self-inductance L on which 
a magnetic flux $e is applied; The total embraced flux $ is related to the 
phase difference between the terminals 4 and 3 

^= ^(<^4-<^3)- (4) 

Four equations (2 ) with the relations (3 ) and (4 ) constitute the coupled 
system of differential equations for the phases ipj and describe the dynamical 
behavior of the 4-terminaI SQUID in the presence of thermal noise. They 

have the integral of motion J2(pj/Rj — C , where the value of the constant 

i 

without loss of generality can be put to zero. Thus we have three independent 
phases and it is convenient to introduce new variables 

1 1 

^2-fl = 0, (^4 - (^3 = 0, -((^1 + (^2) - -{^3 + ^4) = X- (5) 

In the following for simplicity we will consider the symmetric case Ri = 
R2 = R3 = Ri = -R/2. We use the following dimensionless units: current 
in units of Iq = nAl/ (iekTcR), voltage in units of Jq-R and time in units of 
h/{2eRIo). The dynamical equations for the variables 9, 0, x take the form 
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^ = 7 - |sm6' - sin I COS I cosx + Sv0{t), 
# ^ ^ _ 1 sin - sin f cos f cosx + Sv^{t), (6) 
^ = — sin X cos I cos | + Sv-)^{t) , 

where 0e — y^e and JC — ^IqL is the dimensionless self-inductance. The 
"random forces" Sv0,Sv^ , and 5v^ are related to the initial dimensionless 
fluctuating voltages in eq.(2 ) : 



Sv0 = Svi - 5v2, Sv^ = Svs - 5v4, 5v^ = "2 ('^'^i + '^■^2) + 2 ('^'^s + ^'^^a)- (7) 

The corresponding spectral densities in dimensionless units are 

2ekT 

5,, = 5,, = 4r,5,, = 2r,r = ^. (8) 

The system of the coupled equations (6) determines the time dependencies 
of the phases 6{t),(f)(t) and x(^) in the presence of thermal noise at given 
values of the control parameters - the transport current / and the external 
magnetic flux 0e , which are assumed to be independent of time. The time 
derivatives of the phase differences determine the voltages between different 
terminals: 

^21 = 0, 1^43 = 0, \{Vn + ^24) = X- (9) 

3 Current-voltage characteristics in the pres- 
ence of noise. 

We solve the equations (6) numerically by integrating the phases stepwise in 
time with a step At. The Johnson noise voltages 6v{t) are modelled by three 
uncorrelated trains of voltage pulses of constant duration At and amplitudes 
5v. The amplitudes 5v are Gaussian-distributed pseudo-random numbers 
of zero mean value with < Svj >=< Sv"^ >= 2V/At, < 5v'^ >= T/At in 
accordance with the spectral densities (8). The value of F (8) characterizes 
the noise level in our system. 
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In this section we use the above formulated numerical model to calculate 
the noise-affected current- volt age characteristics of the 4-terminal SQUID. 
The observable quantity V, the voltage in the transport channel averaged 
over time, was obtained from the solution of eq.(6) for the phase 6{t) : 

To check the numerical techniques we first calculated numerically the 
noise-rounded IV characteristics for the single resistively shunted Josephson 
junction. In this well known case the phase 6 satisfy 6 = I — sin 6 + 6ve{t) and 
from the corresponding Fokker-Planck equation the average voltage across 
the junction was obtained by Ambegaokar and Halperin |]T5[ analytically. In 
Fig. 2 we compare our numerical results with the Fokker-Planck calculation, 
they are in good agreement. 

Noise-rounded IV characteristics calculated from the full set of equations 
(6) for two values of applied magnetic flux 0e = 0, vr are displayed in Fig. 3 
a,b. We choose L=1.7, F = 0.2, close to the values estimated in experiments 

. Figure 4 illustrates the behavior of current- voltage characteristics in an 
applied magnetic field. For given value of the transport current / and chang- 
ing of the magnetic flux $e the average voltage periodically oscillates with 
amplitude SV . For the detailed analysis of this behaviour in the fluctuation 



free case see ref. \12 



4 Magnetic flux switching in the bistable 
state due to the noise and transport cur- 
rent. 

We now use the numerical model to study the behaviour of the embraced 
magnetic flux in the ring. 

Let us start with the fluctuation free case. For the time independent 
control parameters / and (pe , we can introduce a potential: 

U{9,(j),x\ I,<Pe) = ^ I9-COS ^-cos --2COS-COS-COSX (10) 

such that the dynamical eqs.(6) take the form 
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The stable steady states of the system correspond to the minima of po- 
tential U with respect to variables 9, 0, x given values of / and 0e- The 
minimization of U with respect to x gives that phase x takes the value or 
vr, depending on the equilibrium values of 9 and : cos x = sign{cos | cos |). 
The stable configurations of the noise-free 4-terminal SQUID where studied 
in ref. fl^ . It was shown that at the steady-state domain in the plane of 
(/, 0e) the region exists, inside which the potential (10) has two minima (see 
Fig. 3 in ref.[|l2[). The structure of this bistable state, the height and the 



width of the potential barrier in the 3-dimensional phase space {9, Xi 0) can 
be regulated by control parameters I and 0e- For example, at 0e = tt we 
have two equilibrium values of induced magnetic flux ,which are equally lo- 
cated around vr: 0i,2(-^) = vr± A0. The current circulating in the ring in this 
case equals ±j. The values of the potential U at these two states are equal, 
f/i(/, 0e = vr) = U2{I,4>e = 71"), and in which state the SQUID is located 
depends on the history of the system. It is emphasised that in contrast to 
the case of the usual SQUID, the described bistable state of the 4-terminal 
SQUID exists for each value of inductance £ even for £ < 1. 

The presence of noise, i.e. the three random forces 6v{t) in the dynamical 
eqs.(6), produces small fluctuations of phases 4>,9,x near the equilibrium 
values, as well as transitions between the two states. Each such transition 
corresponds to the switching of the magnetic flux between the two values 0i 
and 02. The switching of magnetic flux must be accompanied by a switching 
of X with a value ±7r or by a 27i slippage of 9. Note, that the phases 9 
and X independent parameters, having different physical meaning. In 
accordance with the eqs.(9) the dynamics of 9 determines the voltage in the 
transport channel, while the time dependence of x means the appearance 
of the voltage between the ring and the transport circuit. In other words, 
the last one means the electrical charging of the ring. The random event of 
one from the two possible processes, switching of ^ or % , is determined by 
the magnitudes of the corresponding potential barriers AUe and Af/^^. They 
depend on the value of the transport current /. By analyzing the shape of 
the potential U{9, 0, x | 0e) surfaces can qualitatively be concluded (see 
ref.|]l3|) that for small values of I the transitions with switching of x ci-re 
more probable than slippage of 9, but with the inreasing of the current / the 
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situation will changes to the opposite. 

To study the processes of switching in the bistable state, produced by 
the thermal noise, we numerically solve the eqs.(6) and obtain the time de- 
pendencies (f>(t),9(t),x{t) ■ Figs. 5 and 6 show some typical traces of the 
phases (f), X and as a function of time t at an external flux value (/)e = vr for 
particular values of the noise parameter F = 0.2 and of the selfinductance 
C = 1.7. 

Fig. 5 shows the behaviour for the transport current 7 = 0. In this case we 
observe the noisy behaviour of phases and the random switching of (p between 
two flux states. When (j) switches, the phase x jumps at the same moment in 
time, with a value of ±7r. In the same time, although the behaviour of 6{t) 
is noisy, no switching is observed in 9. 

In Fig. 6 the transport current equals I — 0.25. It displays the case , in 
which the switching of (p between both flux states may be accompanied by a 
jump of X with a value of ±7r or a 27r phase slippage of 6. In addition, in the 
numerical simulations the case is observed, that x and 6 jump at the same 
time while (p remains constant. 

The behaviour observed in Figs. 5 and 6 corresponds to the above men- 
tioned dependence on the transport current of the potential barriers AUe 
and AU^. The switching of the magnetic flux can be induced not only by 
the noise but also by the transport current /, if it exceeds the critical value 
Ic{4>e)- Fig.7 shows the time dependencies of the phases 4>, x, in the noise- 
free case (F = 0) for the value of the current / = 1.05. This value only very 
slightly exceeds the critical current (see Fig. 3) and we see that the flux 
adiabatically follows ^(t), except the moments of time when the slippage of 
9 takes place and the corresponding switch in the flux occurs. In the same 
time, in this regime, the phase x remains constant. 

5 Summary 

The four terminal SQUID presents a system with three degrees of freedom, 
the phase differences 0, 6, x- In the steady state domain the observable quan- 
tity is the magnetic flux embraced in the ring, $ = ^0. In some region of 
the independent control parameters, the transport current / and the applied 
magnetic flux 0e, the system is in the bistable state with two magnetic flux 
values. The operation of the four terminal SQUID, ie the choosing of one 



8 



from two possible states, can be achieved by variation of / and 0e- At fixed 
values of the control parameters, the magnetic flux switching is produced 
by the thermal noise or by the transport current, which exceeds the critical 
value. 

We have computed the transitions between two flux states stimulated by 
the thermal fluctuations. The switches of the flux are accompanied by the 
jumps in 6' or X and the corresponding voltage impulses. The existence of 
a bistable state in a four terminal SQUID was experimentally confirmed in 
ref.[jl3|. by the direct observations of the magnetic flux behaviour, which is 
similar to the results obtained numerically in Section 4. 
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FIGURE CAPTIONS 

Figure 1. The 4-terminal SQUID. The area of the region closed by 
the dashed lines (the Josephson 4-terminal junction) is of the order of the 
coherence length squared. 

Figure 2. A comparison of the numerically calculated noise-rounded 
IV characteristics (points) for a single resistively shunted junction with the 
analytical results from ref . (solid lines) . The dashed line is the noise- free 
characteristic. 

Figure 3. Current- voltage characteristics of the current-biased transport 
circuit of a 4-terminal SQUID in the presence of noise for F = 0.2, £ = 1.7 
and for two external flux values 0*^ = (a) and = vr (b). The solid lines 
are the noise-free characteristics. 

Figure 4. The changing of the noise-rounded current-voltage charac- 
teristic of a 4-terminal SQUID in an applied magnetic field. SV shows the 
amplidude of the voltage oscillations as a function of the external flux for 
an applied transport current / . 

Figure 5. The three phase differences 0, x s-^^d 6 versus time t for 
(f)'^ = 7T, C = 1.7 in the presence of noise with F = 0.2. The transport current 
/ = 0. Transitions between the two flux states occur (jumps in 0) which are 
accompanied by jumps in x by ±7r. At zero transport current no switching 
in 9 occurs. 

Figure 6. This figure can be compared with Fig. 5 but now the transport 
current is equal to / = 0.25. Again 0^ = vr. This is an intermediate situation 
where the switching of is accompanied by jumps in x by ivr or in 6 by 
27r. At higher currents the x-fluctuation mechanism will eventually be taken 
over by the ^-fluctuation mechanism. 

Figure 7. The time dependencies of the phase differences 0, x and 9 
in the noise-free case (F = 0) for the value of the current / = 1.05, which 
slightly exceeds the critical current and for 0^ = vr. At the moments of time 
when the 27r phase slippage of 9 takes place the switch in the flux occurs. 
The phase x this regime remains constant. 
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